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1. INTRODUCTION 

In accordance with our Technical Monitor this quarterly report combines 
the first two quarters of our LANDSAT-2 research project from the receipt 
of the first data on May 23rd until November 23rd 1975. 

1,1 Participants 


The project is carried out primarely by a joint research group of the 
Department of Geography, University of Zurich (Prof, H. Haefner) and the 
Department of Photography, Swiss Federal Institute of Technology, Zurich 
(prof. W,P, Berg), The following persons are engaged in the project: 


- Klaus Seidel, Ph.D. 

- Robert Binaegger, Ph.D, 
~ Rudolf Gfeller, Ph.D, 

- Arthur Funkhouser, M.S, 

- JUrg Lichtenegger, M.A. 

- Daniel Nuesch, M.A, 

- Rene Muri 

- Karl Stana 

- Hans-Peter Stirnemann 

- Fritz Fasler, M.A, 

- Anton Paschke 

- Otto Matt 


project leade] 


system specialist 
system specialist 
system specialist 


(part time) 

(until March 31st, 1975) 

(since November 1st, 1975) 
(part time) 

(part time) 

(part time) 

(part time) 

(part time) 

(part time) 

(part time) 


Dr. Klaus I tten, having spent one year as a visiting researcher at GSFC in 
1974 (ESA‘.*grant) returned to the Department of Geography, University of 
Zurich, and is continuing his specific research on the comparison of 
different existing image processing systems to operationally classify 
snow covered areas (3.l). 

In addition the following institutions participated in the research project; 

- Department of Geography, University of Benne (Prof ,B,Messerli/Ds',,MiUiniger) 

- Swiss Institute of Snow and Avalanche Research, Veiasfluhjoch (Prof. iM.de Quervain) 

- Swiss Meteorological Institute (Dr. A. Piaget) 
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1,2 International contacts 


Various international contacts could be establisJied or strengthened during 
the reporting period t 

Dr, K. Seidel participated in the "4th Annual Remote Sensing of Earth Resources 
Conference at the Space Institute of the University of Tennessee/'I'-ullahoma" 
from March 24th to 26th, 1975 and visited several other institutions in the 
US, such as Office of Remote Sensing, Pennsylvania State University; EROS 
Applications Assistance Facility, Bay St. Louis, Miss,; Jet Propulsion 
Laboratories, Pasadena, Calif, etc. 

Prof. H. Haefner attended the "Earth Resources Survey Symposium" from June 8th 
to 12th, 1975 in Houston and visited the Space Technology Center, University 
of Kansas, Lawrence and the L.B. Johnson Space Center in Houston. 

Dr, K. Itten presented the findings of his research (5.l) at the "Workshop 
on Operational Applications of Satellite Snowcover Observations" from August 
18th to 20th in South Lake Tahoe, Calif, and visited the Stanford Research 
Institute (Dr. W. Evans) , the School of Earth Sciences, Stanford University 
(prof .R, Lyon) and the GSPC. 

The Department of Photography, Svriss Federal Institute of Technology, Zurich, 
organised a Symposium on "Digital Image Analysis, which attracted quite a many 
researchers not only from Switzerland but from abroad as well. 

At the "10th International Symposivun on Remote Sensing of Environment" in Ann 
Arbor from October 6th to lOth, 1975 Mr. D. Nuesch presented a status report 
by Prof. H. Haefner et al. on the applications of remote sensing in Switzerland 
with special reference to our LANDSAT -activities. 

Permanent contacts and exchanges of ideas and results regarding our LANDSAT- 
studies are taking place with the "Euratom Ispra" (Dr, B, Sturm) and the 
"Deutsche Forschungs- und Versubhsanstalt fUr Luft- und Raumfahrt", Oberpfaffen- 
hofen, Germany, 


2 . TECHNIQUES 
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2,1 Data coverage 


For a complete coverage of Switzerland two frames each on three following 
days are necessary (Fig. l). During the reporting period the material listed 
in Fig. 2 was received. 

Prom this tabulation it is easy to recognize that the coverage is very iidiomo- 
geneous and not well suited for studies of dynamic processes such as changes 
in snow cover / land use etc. This fact vras brought to the attention of NASA 
by letter of September 6th, suggesting that 

- both frames of a daily cycle should be taken together always, 

- preference should be given to the easterly orbit if possible (to reach a 
more continuous coverage of our main test areas in Grisons), 

- a complete coverage of Switzerland (on three following days) would be of 
great interest (since never achieved with LAND SAT-1 and 2 until now). 

2.2 Methodological aspects 


Digital as well as analog processing of the LANDSAT-images was undertaken, 
but with a clear priority to the digital methodology, because of their greater 
flexibility. 

The method including all preprocessing steps for geometric and radiometric 
corrections as well as for the output on film as developed during our LANDSAT-1 
project was improved further (3.l). It may be summarized as follows: 

To classify the data a multivariate procedure, the linear discriminant analysis 
was used. Prom well considered training groups criteria for the identification 
of the sites were computed stepwise, with given level of tolerance and each 
picture elements assigned to just one of the training groups. The results were 
displayed on "single-symbol-separation-films" used for color overlay production. 

For snow mapping all four MSS-channels are used as variables whilst for classi- 
fying land use features (3. 2 ) six synthetic ratio variables were used in 
addition. 
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Tho accuracy of our method depends prifflaroly on the selection of representative 
training c-elis. Therefore extensive ground observation and measurements were 
carried out (2. 4/2. 5) for our snow mapping studies. The ultimate goal would 
be to establish a comorehensive collection of samples, to minimize or even 
abolish ground truth in the future. 

Fig. 5 summarizes the classification system and the computer programs for 
digital picture processing of LANDSAT-data as available today. 


2.3 Aerial underflights 


During several LANDSAT-passes underflights could be arranged. Unfortunately 
the satellite was not always covering our areas when a flight was undertaken. 
Therefore only for two dates matching images of LAHDSAT and high altitude 
aircraft could be achieved. The technical data are the following: 

- Camera: Omera . 

- Objectives; f = 44 mm and 100 mm 

Film; panchromatic 

- Quality of photography: good 

- Altitude: approx, 12 '000 m above sea level 

- Scale: approx. 1:225 '000 and 1:95 '000 

- Dates: April 22nd and August 8th, 1975 

2,4 Ground truth 


For all LAl^DSAT passes with good weather conditions extensive ground surveys 
were undertaken in the main test area of Grisons (between Davos and Lenzer— 
heide), surveying and registering the following objectives; 

- observation of position of temporary snow line (visible observations, 
photographs, mapping) 

- observation of various snow-parameters at selected sampling points: 

= type of snow 

== characterization of snow surface (forms, cleanness. Wetness, age etc.) 

= characterization of snow pack (granularity, wetness etc.) 
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= snow temperature (at surface) 

= snow depth 

= snow mass 

== air tempei'ature 

- air moisture 

= slope angle and exposure 

= cloudiness 

= daytime 

This information will le of great assistance in selecting the "training cells" 
for the digital snow classification. 

2.5 Measurement of spectral reflectance of snow 


In addition an extensive measuring program was set up to get more information 
on the spectral nroperties of the snovf surface. An EXOTMCH-100 radiometer was 
used for this purpose on the ground as well as from a helicopter. During the 
melting period 1975 (February to June) measurements were carried out for 
about 20 different days on the ground and for two days from a helicopter in 
the vicinity of Davos. 

Measurements under various conditions (sun angle,, slope angle, exposure, sun, 
shadow etc.) were undertaken for six different snow types. In addition spe- 
cific measurements for "dirty" snow covered partly with dust and sand were 
taken. At the same time ground observations as described in 2.4 were carried 
out at the measuring points. 

Again the data will be processed further to improve our classification 
algorithms. 


2,6 Laboratory equipment 
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For our digital processing the following specific equipment was at 
disposition: 

— computer 

== CPC 6400/6500 with standard peripherals 
= IBM 370 vrith standard peripherals 

- data output 

= OPTRONICS PHOTOMION SYSTEM P I7OO 
= computer prints 


3 . ACCOMPLISHMENTS 


3*1 Snow napping 

Methods and results of our snow mapping project were reported in the BRTS-1 
final report. 

Regarding snow mapping emphasis was layed during the reporting period on 
_ evaluation of EREP-S-192 digital data for separation of snow and clouds 
improvement of methodology (2,2) 

- careful selection of training cells by ground observation (2.4), ground 
measurements (2.5) and aerial under f lights (2.3) • 

It is intended to continue snow mapping vrith the improved classification system 
to investigate changes in the areal extent of the snov7 cover. But for this 
purpose additional data of the same test area (2,l) is essential.' 

Interesting results occui-red from the research by Dr, }C, Itten on the feasi- 
bility of three different existing image processing systems to operationally 
classify snow coVei'ed areas* 


1) HABPNER^ H. : Snow Survey and Vegetation Growth in High Mountains (Swiss Alps) 
and Additional ERTS-Investigations in Switserland, Type II Final Report, 
January 1975. 
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LARSYS Version 3 was used first to obtain a digital snow cover classification. 
V/ith the supervised, statistically proven method, seven major ground cover 
classes vrore distinguished. 

The overall achieved classification accuracy regarding to the training areas 
amounted to 92 fo. The process with LARSYS Version 3 is stepvrise, well defined, 
and tho res'alt is - good ground truth given -* q,uito accurate. But the whole 
procedure is very time consuming and in that sense costly. It is therefore 
very questionable whether to use this approach for operational snow mapping 
at all. 

Having gained much experience with the data, the STANSORT-2 system of the 
Stanford University was used to check the feasibility of a highly interactive 
semi-supervised clustering as an approach to the snow mapping problem. V/ithin 
very short time, by comparing intermediate results with different density 
slicing and cluster tolerance factors, v/ith known test area patterns, a 
classification very similar to the one obtained with LARSYS could be achieved. 
The fast evaluation of intermediate results vie the b&w TV-screen was a main 
advantage and time saving factor. The simple and thus fast cluster algorithm 
proved to be very successful in this specific application. 

The third and probably one of the most advanced systems in operational image 
processing tested, was the Image-100 of General Electric. The same COT' s and 
test area, and part of the same ground truth were used again. A classification 
into the same classes v/as achieved in an extremely short time. The simple and 
fast parallel epiped classifier was used, and the classv/ise intermediate 
results were controlled immediately on the color TV-screen. With Image-100 
only 1 or 2 training areas v^ere used for each ground cover class. This to- 
gether v/ith the relatively coarse classifier is an explanation for, the some- 
what different results compared to the classifications by LARSYS and STA.HSORT. 

The results of the comparative study are summarised in Pig. 4. 


3.2 Land use mapping 


Significant rosulta could 'bo gained fi*om the research on land use mapping 
with the general aim of applying LALT)SAT“MSS-data as information source 
to thematic mapping. 

The test site is an area of approx. 26 x 39 1cm in the Po-Valloy near Milano. 
It is almost flat land which means that no influences by relief are present 
and not considered in this study. The overpass of October 7th, 1972 (B-1076-^ 
09442) provided again the best suited data. 

The following land use categories were considered: 


Land Use category 

(1) V/ater 

(2) Built-up areas 

(3) Fields, cultivated with dead 
vegetation and bare soil 

(4) Fields, cultivated vdth living 
vegetation 

(5) Forestry 


abrevation for group 

(as listed in figures 5-0-7) 

SEE 

GRO, KLE, HEL 
TOT 

LBB 

¥AL 


A careful selection and delineation of the training cells vras achieved by 
field observation as well as statistical tests to reach homogeneous groups 
representing a single land use category. To characterize a special category 
it was sometimes necessary to use different subsets (e.g. for built-up 
areas). 

Applying a stepwise discrimination analysis (BJID07M) it was possible to 
weight the variables and to determine the most suitable ones for a separation 
of the different categories, which is used first. A classification based 
strictly on the 4 original KSS channels as variables showed that only channel 5 
and 7 v/ere needed for an optimal separation (Fig. 3 ) » whilst by adding 6 ratio 
variables the following succession in the use and importance of the variables 
was achieved; 7/5 - 6 - 7/6 - 4 (Fig. 6 ). 
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Tho results show that "by using the original LAIfD SAT -variable o only the 
categories "TOT" and built-up areas ((GRO, KL3^ II3L) were separated with 
the lowest accuracy, but that by adding the ratio variables an improvement 
of about 7 could be achieved. On the other hand the separation of the 
built-up areas (GHO, HSL, KLE) against vmter (SEE) worsened. 

Therefore the use of ratio variables does not always improve the classification, 
which consequently cannot be done in a single run. It is always necessary to 
experiment v;ith slightly changed initial choices, both for variables and 
learning sets used. According to each task, appropriate choice will ensure 
optimum sharpening of the classification groups and distinct separation 
between them. 

Having established the training cells and the* algorithms the entire test area 
was classified. Fig. 7 shows the corresponding classification matrix. 

Finally the results vrere reproduced in cartographic form (scale 1:200*000) 
with the photomation film -write system. By preparing color overlays (instead 
of black and vfhite prints) a clear presentation of the different categories 
could be gained. For this purpose three black and white overlays were printed 
for: 

1) Categories GRO/SEE/WAL cyan 

2) Categories GRO/KLE/HEL — — magenta 

5) Categories KLB/WAL/LEB — yellow 

By combining the three overlays, a color composite was prepared (Fig. s) 
showing the different land use categories in the color specified in Fig. 7 , 

Fig. 8 shows the final classification Of the test site and the output in 
form of a colored thematic map. 

The results were verified in the field and vrith existing topographic and 
land use maps. 
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4. SIGITIPIOARr iiSSULTS 


Sii;nificant results could be aerxieved in the rnethodology of snow mapping 
and in land use mapping from digital KSS-data. 

Kethodolo^?/ of snow manning 

Applying tho same LAlTDSAT-data to three different image processing systems 
(LARSYJ, Vers. 3; STANSORT-2 and GE-IMAGS~100) the following results for 
snow mapping in forested mountainous terrain could be gained: 

Generally observable is a tradeoff between the classification accuracy and 
the time or money needed in the three different approaches. VJhile the 
control and accuracy achieved with LAHSYS is remarkable, time and effort 
to perform the processing favourize the systems S'TAl’J'SORT and IMAGE-100. 

A substantial increase in the applicability of satellite data to operational 
snow mapping could be demonstrated by connecting a multispectral approach 
vrith digital image processing techniques. It is believed that for the future 
a highly interactive, specially designed system together with a skilled 
applications expert vail provide maximum operational use in satellite snov/ 
cover observations. 

Land use maonina: 

A system was developed to operationally map and measure the areal extent of 
various land use categories for updating existing and producing new and actual 
thematic maps shov?ing the latest state of rural and urban landscapes and its 
changes. 

The processing system includes: 

- preprocessing steps for radiometric and geometric cori'ections 

- classification of the data by a multivariate procedure, using a stepvdse 
linear discriminant analysis (BICD07M) based on carefully selected training 
cells. In addition to the 4 MSS-bands synthetic ratio variables are used 

to improve the classification. The importance of the variables (in accordance 
to their order of succession considered in the classification)was ; 

7/5 - 6 - 7/6 - 4. 
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- output in form of color jaeps by printing black and white thomo ovorlaya 
of a cclccted scale (e.g. 1:200 ’OOO) with an Optronic.s photomation system 
and its coloring and combination into a color composite. 

'fhe method was applied to a test area in the Po-Valley near Kilano, Italy. 
5. PUBLICATIONS 


The following publications were already forvrarded to NASA: 

1. BINZ3GGSR, R. : ERTS-I-Iultispektraldaten als Informationsquelle fiir 
thematische Kartierungen, Diss. Univ. Zurich, 1975- 

2. GFSLLER, H. : Untersuchungen sur automatisierten Schneeflachenbestimmung 
mit Kultispektral-Aufnahmen dos Erderkundungssatelliten ERTS-1, Diss. 

Univ. Zurich, 1975* 

5. HASFInSR, H. + JiBSSERLI, B. : Erderkundung aus dem ¥eltraum - Das schwei^ 
serischen ERTS- und EHEP~Satellitenpro jekt, Geographica Helvetica 3/75. 

4. WINIGER, M. ; Untersuchung der Nebeldecke mit Hilfe von ERTS-i-Bildern. 
Geographica Helvetica 3/75 • 

5. GFELLBH, R. + SEIDEL, K. : Automat isierto Schneeflachenbestimmung mit 
digitalen Multispektraldaten des ERTS-1, Geographica Helvetica 3/75. 

6. HAEFNER, H. + GEISER, U. ; Kartierung von Hbhengrenaen zwischen Ht. Blanc 
und Gotthard-Massiv mit Skylab-EREP-Aufnahmen, Geographica Helvetica 3/75. 

Publications in print; 

1. SEIDEL, H. + GFELLER, R. + BINZEGGER, R. ; Snow and Vegetation Classification 
by Means of Digital LANDSAT-M3S-Data, Proc. 4th Ann, Remote Sensing of Earth 
Resources Conf,, Univ. of Tennesse, Tullahoma (l975). 

2. ITTEH, K. : Approaches to Digital Snow Mapping with LANDSAT-1 Data, Proc. 
V/orkshop on Operational Applic. of Sat. Snowcover Obs*, South Lake Tahoe (l975) 

3. HAEPI^ER, H. + ITTEN, K. + SCHALDA, E. ■+ LOTGER, M. + SEIDEL, K. ; Applications 
of Remote Sensing in Switssrland, Proc. 10th Int. Remote Sensing of Env. 
Symposium, Ann Arbor (l975). 


- 12 


6. PROBLEMS AN!) RECOI>mEI®ATIONS 


The main problem for our project is as already mentioned (2.l) the un- 
systematic and insufficient coverage of our test site, which does not 
permit to map the actual changes of the extent of the snovr cover during 
a certain time period. The necessary recommendations as stated in 2,1 
are to give priority to the easterly orbit of Switzerland (Fig. l) and 
to always take data of both scenes during a single orbit. 
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Fig. 2: LANDSAT-2 scenes received during reporting period (index as in Fig, l) 
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Fig» 3s Computer programs for digital image 
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-MSH - COMPUTATION. OF SCANLINEFACfORS 
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4: Comparison of classification results (number represent areas 

in percent), (after ITTEN, 1975) 


Cover Type LARSYS Ver. STANS 0!RTi^2 ; GE Imager-100 


Dry Snow 
Metamorph . Snow 
Forest. w.Snow 
Interzone 
Bare Porest/Veg, 
Shadow -f Water 

31.8 

22.1 

27.2 

9.4 

6.2 

0.4 

30.8 
22.5 

27.8 
11.2 

6.0 

0.2 

30.9 

21.1 

27.1 

4.6 

10.4 

0.2 

Total Snow 
Covered Area 

85.8 

86.7 

81.4 

Total Area 
Bare of Snow 

11.3 

11.8 

12.9 

Unclassified Area 

2.9 

1.5 

5.7 

Accuracy/Test fields 

92 (calc. ) 

90 (est. ) 

87 (est. i 
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5 : Separation of training cells using the 4 MSS-hands as variables 

(in frames are subgroups, which together form the main land use 
categories), (after BINZSGGER, 1975) 


step 1 
(variable 7) 
category 
pixel 

total pixels classified 
of each category 
SEE GRO KLE HEL TOT LEB WAL 

% of correct 
classification of 
each category 

SEE 

GRO 1’067 
KLE 1*025 
HEL, 117 
TOT 998 
liEB 1 ' 002 
WAL 736 

rill2l 4 0 "o1 0 O 0 

96,5 

IToo,^ 
1 97,01 

I 88,0 

1 92, ol 

{ 30 993 44 0 

1 O 257 631 15 

1 O 1 8 SO 

0 0 0 
122 0 0 

97.0 

88.0 
59,0 

O O 113 412 
O 0 O 78 
0 0 25 393 

458 15 ^ 

3 921 ol 

47,5 

92,0 

306 12 1 0l| 

0^0- 


step 2 


(variables 
7 and 5) 
category-*- 
:pixel 

total pixels classified 
of each category 

SEE GRO KLE HEL TOT LEB WAL 

^ of correct 
classification 
each category 

SEE 116 
GRO 1 • 0c7 
KLE ; 1*025 
HEL 117 
TOT 998 
LEB 1*002 
WAL 736 


97,5 

fIoo,o] 

100,0 

87,5 

r84,T^ 

,100,0 

' 1 947 119 0 1 O O 0 

1 . 0 186 713 0121 1 4 

1 0 0 10 97 1 lO 0 0 

IQO, o 

87.5 

91.5 

0 0 144 8 [1777 25 44] 

Q 0 0 o| 11 937 541 

80,5 
94, 5 


96,5 

MIV 1*011 

FLU 341 


80,5 
41, o' 

I3S13 

&il! 










Fig. 6; Separation of training cells using the 4 MSS-hands and 
synthetic ratio variables, (after BINZEGG-ER, 1975) 


step 1 
(variables 
7 / 5 ) cate- 
go ry+ pixel 


total pixels classified 
of each cateogry 

SEE GRO KLB HEL TOT LEB WAL 


SEE 

116 

I|l02| 

12 

2 

0 

0 

0 

0 

88,0 

GRO 

1'067 

115 

650 

62 

238 

2 

0 

0 

89,0 

KLE 

1*025 

0 

97 

504 

326 

95 

0 

3 

90, 5 

HEL 

117 


46 

16 

52 

_ J3 0 

0 

7^ 

97,5 

TOT 

998 

0 



"”3T 

620 

2 

62,5 

LEB 

1*002 

0 

0 

0 

0 

7 

799 

196 1 

80,5 

WAL 

736 

0 

9 

0 

0 

41 

95 

[6^ 

81.5 


% of correct 
classification of 
each categoi'y 


100,0 

90,51 

I 69,5, 

100,0 


sxep 

(variables 

total pixels classified 

^0 of 

correct 

7/5-6-7/6-4) 

category •V- 
pixel 

of each category 

SEE GRO KLE HEL TOT LEB 

WAL 

classifxcation 
each category 

SEE 

116 


11 0 0 

0 0 

0 

90,5 


100,0] 

GRO 

1*067 

21 
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Fig. 7: Classification matrix for land use mapping of total test are/i 

(including color of category represented in Fig. s). 

(after SEIDEL/GPELLER/BINZEGGER, 1975) 
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